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Abstract. It is shown that taking into account a recoil for radiation scattering in the Ly- 
a line can lead to a noticable acceleration of primordial hydrogen recombination. Thus for 
ACDM model a decrease of ionization degree exceeds 1% for redshifts z in a range 800 - 1050 
achieving ^ 1.3% at 2; = 900. Corresponding corrections to the cosmic microwave background 
power spectra can achieve 1.1% for TT spectra and 1.7% for EE ones. Radiative transfer in 
these calculations was treated in a quasistationary approximation. Numerical solutions are also 
obtained in diffusion approximation for a nonstationary problem of Ly-a line radiative transfer 
under partial frequency redistribution with a recoil. An evolution of a local line profile is traced 
to as well as an evolution of a relative number of uncompensated transitions from 2p state down 
to Is one. It is shown that taking into account nonstationarity of Ly-a line radiative transfer 
can lead to an additional acceleration of primordial hydrogen recombination. 
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INTRODUCTION 

In connexion with the proposed considerable increase of CMBR observations precision the 
requirements grew as well to predictions accuracy of the theory being used to interpret the 
observations. In particular this concerns the theory of cosmological recombination. Required 
accuracy is about 0.1%. So at last time a few papers appear (Dubrovich, Grachev, 2005; 
Chluba, Sunyaev, 2006; Kholupenko, Ivanchik, 2006; Rubino-Martin et al., 2006; Chluba et al., 
2007; Wong, Scott, 2007; Chluba, Sunyaev, 2007) where recombination history of primordial 
hydrogen is calculated with a proper allowance for some fine effects which however are essential 
at the required level of accuracy. 

In this paper we study an effect of allowance for a recoil under radiation scattering in 
Ly-a line on recombination history of primordial hydrogen. Earlier we already carried out 
calculations of such kind (Grachev, Dubrovich, 1991) for purely baryonic models. An effect of a 
recoil turned out to be practically negligible. However our calculations fulfiled in the framework 
of present-day ACDM models show that for these models characterized by low baryon density 
the effect of recoil becomes nearly ten times larger. The reason is in the fact that at low 
baryon density the relative role of Ly-a photons escape from the process of scattering due 
to Universe expansion (so called "intrinsic escape") considerably increases as compared with 
another main mechanism of irreversible recombination, namely, two-photon transitions from 
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the second level. Therefore the role of fine effects in treating of radiation transfer in the Ly-a 
line grows correspondingly. 

The second curcumstance to which we pay attention is an influence of nonstationarity 
of radiaton transfer in the main resonance line on recombination history. Up to now a quasi- 
stationary approximation is in usage. Namely, when computing hydrogen recombination history 
one makes use of stationary Sobolev approximation for the number of uncompensated transi- 
tions from 2p state down to Is one (so called net radiative bracket (NRB)). However it would 
be very important to determine characteristic time necessary for NRB to reach a steady-state 
limit. In this connexion we fulfiled model numerical calculations of time dependent radiative 
transfer in a resonance line under partial frequency redistribution (PFR) allowing for a recoil 
in the process of scattering. It turned out that time of reaching a steady-state limit for NRB 
may well exceed 1% (in Az/z units) in the range of redshifts near z ~ 1000. 

RECOMBINATION DYNAMICS ASSUMING QUASI-STATIONARY RADIATIVE 

TRANSFER 

It is known that Sobolev approximation widely used in astrophysics to treat radiative trans- 
fer in spectral lines in moving media gives an exact solution of a stationary problem for the 
case of two-lewel atoms in an infinite uniform expanding medium assuming complete frequency 
redistribution (CFR) in a single scattering. Then for a relative number of uncompensated 
transitions (so called NRB) we have 

r2i = [n2(A2i + S21 J12) - n,BuJi2]/n2A2i = /5(1 - B/S), (1) 

where 

/3 = 7 (1 - e-i/^) (2) 
has a sense of photon escape probability out of the process of scattering, 

^ A?2^l[(c/2M)^l -^2] 

is a dimensionless velocity gradient proportional to velocity gradient H (which is the Hubble 
factor in cosmological context). A quantity rs = I/7 is called Sobolev optical distance. It is 
an optical depth of an infinite expanding medium along photon path. Above we use standard 
designations for Einstein coefficients of transition probabilities and for statistical wheights of 
levels. Moreover J12 is the mean intensity of radiation in a line, A12 is a transition wavelength, 
S = {2hc/ Xi2)[gin2/ {g2ni — ^'1^2)] is a line source function, B is the Planck function at the line 
center (we suggest that far in the violet wing of the line radiation intensity is planckian and 
we also ignore frequency dependence of the Planck function within the line). From eqs. ([2]) 
and ([3]) we have equation f3 oc H/Qb which clearly demonstrates a grouth of escape probability 
when f^B decreases. 

As it was shown by Chugaj (1980) rejection of CFR i.e. usage of PFR does not affect on 
a mean number of scatterings and on excitation degree of atoms. Therefore if one neglects 
induced radiation and suggests that 7 -C 1 then the equations given above may be used for 
PFR (Grachev, 1989; Grachev, Dubrovich, 1991) if only recoil under scattering is neglected. 
Allowing for recoil leads to the growth of photon escape probability due to additional redshift 
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of photons. From analytical solution of the problem in diffusion approximation (Grachev, 1989) 
the following expression is obtained for it suggesting zero boundary condition in a far violet 
wing: 

^ = z(-oo)7, 7<1, (4) 
where correction factor i{—oo) > 1. It depends on parameters 

p = 26x^, a = x^/xx, (5) 

where 5 is a recoil parameter and x^ and xx are characteristic dimensionless frequencies: 

S = hiyi2/Mvti,c, x^ = (3aA/27r7)i/^ xx = y^A/2(l - A). (6) 



1 -1 

, a < 0.38, (8) 



Here M is an atom mass, Vth = y 2/cTe/M is a most probable velocity of atoms thermal motion, a 
is the Voigt parameter, A is an albedo per single scattering (a probability of photon "survival" in 
a single scattering). Characteristic frequency x^ was introduced by Chugaj (1980). It separates 
domain of photons diffusion (|x| < x^) from domain of drift (|x| > x^) due to differential motion 
(expansion) of the medium (here and so on x = — z/i2)/Az/d, Az/d = i^uVth/c). Results of 
numerical calculations of i{—oo) can be approximated by equations 

i(-oo) ^l+p/[l + (aV3)(4-a2)/(2 + (T2) + frV6], a > 0.38, (7) 

and 

i(-oo) ^ 3 / e~^y~y y^dy 
. Jo 

with an error less than 0.3%. It should be noted that for i{—oo) — 1 an error of these equations 
is less than 8%. 

For non-zero boundary condition /(+oo) = B in a. far violet wing eq. ([1]) turns then into 
(Grachev, Dubrovich, 1991) 

r2i=i[i{-oo)-B/S]. (9) 

Computations of hydrogen recombination history using eqs. ([7]) - (jH]) for 2p - Is transition 
were carried out by us earlier in the framework of 60-level model of hydrogen atom for purely 
baryonic universe models and in particular for values of parameters being used in a pioneer 
works on cosmological recombination (Zeldovich et al., 1968; Peebles, 1968): f2tot = 1, = 1, 
Hq = 98 km/s/Mpc, Tq = 2.7 K, Ine = 0. Results of these our calculations are shown in Fig. 
1 along with the results for ACDM model with parameters Qtot = 1) = 0.7, I^dm = 0.26, 
= 0.04, Ho = 70 km/s/Mpc, Tq = 2.728 K, Fhc = 0.24. Namely, in Fig. 1 there are 
graphs of relative variation (in percents) of a fractional electron density Xe = Ue/n-a due to 
allowance for a recoil under Ly-a photons scattering. (Here rie and tt-h are the number densities 
of electrons and of hydrogen atoms and ions respectively.) As a mechanism of Lyman photons 
"destruction" per scattering we consider excitation and ionization of atoms from 2p state due 
to absorption of background black-body radiation. So we have for albedo per a single scattering 
an expression 

A = ^2p,ls ^ ^^2p,ls + R2p,c + 51 -^2p,i^ • (10) 
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Figure 1: Variation of hydrogen recombination history due to effect of a recoil under scattering 
in Ly-a hne. Parameters values: 1) fitot = 1, = 0.7, fie = 0.04, fioM = 0.26, Hq = 70 
km/s/Mpc, To = 2.728 K, Y^e = 0.24; 2) fi^ot = 1, ^^a = 0, = 1, ^^dm = 0, i^o = 98 
km/s/Mpc, To = 2.7 K, Yne = 0. 



One can see in Fig. 1 that for a present-day Universe model recoil contribution increases more 
than ten times and is in the level of 1% near the surface of last scattering. 

Variation of primordial hydrogen recombination history has an influence on the power spec- 
tra of CMBR angular fluctuations in intensity (temperature) and polarization. Corresponding 
graphs are shown in Fig. 2. Computations were fulfiled by means of CMBFAST package (Sel- 
jak, Zaldarriaga, 1996). It is seen that relative variations can reach 1.1% for TT spectra and 
1.7% for EE spectra for multipoles / < 2500. 

NONSTATIONARY RADIATIVE TRANSFER in Ly-a LINE 

To clear up how fast NRB reachs a steady-state limit we have solved a model problem of 
nonstationary radiative transfer in isolated spectral line. It is suggested that at some initial 
moment of time [t = 0) there is no radiation and at this moment primary sources of line 
radiation with a constant power are switched on. The main integro-differential equation of 
radiative transfer has a form 

dlit^ _ ^ dl{t,x) ^ j^^^ ^) ^ ^ /•+°" ^(^^ ^y^f ^ _ A)0(a;)5o, (11) 

Ot ox J-OD 

where I{t, x) is a radiation intensity {x = {v—vi2)//S.vd is a distance from a line center measured 
in Doppler widths; t is a time measured in units of mean time between succesive scatterings), 
R{x, x') is a frequency redistribution function per scattering, (f){x) is an absorption coefficient 
profile, 7 is a dimensionless velocity gradient (see above eq. ([3]) in which however a contribution 
of induced radiation should be neglected), A is an albedo per a single scattering. The following 
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Figure 2: Variation of CMBR power spectra of temperature (TT) and polarization (EE) fluc- 
tuations. Computations are fulfiled with the CMBFAST package (Seljak, Zaldarriaga, 1996). 
Values of parameters are the same as for the curve 1 in Fig. 1. 

normalizations hold 

/+00 r+oo 

R{x,x')dx = (j){x'), / (j){x)dx = l. (12) 

-OO J — oo 

Initial and boundary conditions for eq. ffTT]) are /(0,a;) = 0, /(r, +oo) = 0. 

In the context of cosmological recombination an absorption coefficient profile is the Voigt 
one with parameter a = y42p,is/47rAi/D. As concerned frequency redistribution function with 
allowance for a recoil it was derived in the paper by Basko (1981) where diffusion approximation 
was also obtained for the integral term (in transfer equation) describing an act of a single 
scattering. Using this approximation one can reduce eq. (ITT]) (with 5*0 = 1) to the form 

dl dl / 2 2p\ I(T,y) If 2\dl 1 d^I a\ ^, , , , 

where new dimensionless frequency and time are y = x/x^, t = ^tjx^ and e = a/nx^ is a small 
parameter. Parameters p and a are defined above by eqs. ([5]) and ([6]). When obtaining eq. 
f|T3l) the profile of absorption coefficient was taken approximately in the form (Grachev, 1989) 

(j){x) = 5{x) + a/nx'^, (14) 

where the first term (J-function) describes Doppler core and the second one corresponds to the 
Lorentz wings of the Voigt profile. 

In a steady-state case (when dl/dr = 0) analytical solution of eq. f[T^ was constructed by 
Grachev (1989) and then has been used (Grachev, Dubrovich, 1991) to compute recombination 
history of primordial hydrogen (see also above in the preceeding section). In the present paper 
we obtain numerical solutions of nonstationary eq. (fT3|) for the values of parameters typical 
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Figure 3: Evolution of a local Ly-a line profile in the case of non-conservative scattering. An 
upper profile corresponds to r = cxd. Values of model parameters are the same as for the curve 
1 in Fig. 1. 

for cosmological recombination. (As a matter of fact the function f{r,y) = I{T,y')dy' was 
found at first from equation which is easily obtained by integrating eq. f|T3l) .) We make use 
of numerical method of solving nonstationary problems of radiative transfer theory proposed 
by Grachev (2001). Results are in Fig. 3. The values of ACDM model parameters are the 
same as for the curve 1 in Fig. 1. As for parameters depending on redshift z their values were 
calculated by means of our code which computes hydrogen recombination history (Grachev, 
Dubrovich, 1991). The point z = 1180 near the surface of last scattering was taken. At this 
point we have = 0.300, 1 - A = 2.00 ■ 10"^ a = 8.30 ■ 10-^ 7 = 2.26 ■ 10"^ = 56.0, 
e = 5 ■ 10~^, p = 0.0500, a = 1.12, i(— 00) — 1 = 0.0296. For comparison in Fig. 4 results 
of similar calculations are shown but for cr = which corresponds to conservative scattering 
(A = 1) (then i{—oo) — 1 = 0.0455). (Since in this case eq. f|T3|) becomes homogeneous we 
normalize its solution to unity at r = 00 and x/x^ = —00.) One can see that in both cases the 
profiles approach asymptotically to a steady state limit which corresponds to the upper curves 
on Figs. 3 and 4. 

For NRB one can obtain directly from eq. flTT]) the following expression 
where 

/+00 
(f){x)I{T,x)dx. (16) 
-00 

According to eq. (fT4l) it follows from here that J(r) = /(r, 0) + 0{e). Fig. 5 displays that 
approaching of r2i to the stationary limit (dashed horizontal lines) is going on quite slowly 
especially for non-conservative scattering {a = 1.12). The point is that for non- conservative 
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Figure 4: The same as in the preceeding Fig. but for conservative scattering. Numbers near 
the curves are the moments of time r. 

scattering the hne profile has extended quasi-lorentzian wings (one can show that /(r, y) ~ 
(T^T[3y(y + t)]'^ for \y\ ^ 1, |y + t| » 1) which evolve very slowly towards the steady state 
limit (J(oo,t/) ~ a"^ /2>y) for y ^ 1 (see Fig. 3). According to Fig. 5 the time to reach steady- 
state limit can be as large as r ^ 10 and even more. An interval of redshifts A;z corresponding 
to this time can be estimated from equation dz/dt = —(1 + z)H{z) from which it follows that 



From here one can obtain that Az/z — lA • 10~^At for given above values of parameters at 
z = 1180. Thus non-stationarity of radiative transfer in Ly-a line can lead to an appreciable 
acceleration of irreversible recombination of primordial hydrogen. 

From physical point of view an increase of NRB due to non-steady state radiative transfer 
takes place because the line does not reach a saturation in some frequency regions after primary 
sources switch on, i.e. the are no enough resonance photons in these regions to compensate in 
a maximum degree radiative transitions down by those of up. 



In this paper we study an effect of a more accurate treating of radiative transfer in Ly-a 
line on a calculated recombination history of primordial hydrogen. An urgent necessity of such 
kind investigation is connected with increased role of Lyman photons escape from a scattering 
process in the framework of present day ACDM models characterized by low baryon density. 
We show that alowance for a recoil under scattering in Ly-a line can lead to an appreciable 
acceleration of primordial hydrogen acceleration. As a main one we consider a model with 
the following values of parameters Qtot = 1, = 0.7, Qdm — 0.26, Qb — 0.04, Hq — 70 




(17) 



CONCLUSION 
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Figure 5: Ly-a photons escape probability (or NRB) normalized to Sobolev value. Dashed 
horizontal lines correspond to the steady-state limit allowing for a recoil (upper line: cr = 0, 
lower line: a = 1.12). 



km/s/Mpc, To = 2.728 K, Yne = 0.24. For this model a decrease of hydrogen ionization degree 
exceeds 1% in the range of redshifts z = 800 — 1050 reaching ^ 1.3% at 2; = 900. Corresponding 
variations in computed power spectra of CMBR angular fluctuations reach 1.1% for TT spectra 
and 1.7% for EE spectra for multipoles / < 2500. It should be stressed that radiative transfer 
in these calculations is treated (as well as in other works on cosmological recombination of 
hydrogen) in a quasi-stationary approximation. In order to study an applicability of this 
approximation we have obtained numerical solutions of nonstationary problem of radiative 
transfer in Ly-a line assuming partial (true) frequency redistribution with a proper allowance 
for a recoil. An evolution of a local line profile is traced through as well as an evolution of 
number of uncompensated transitions from 2p state down to Is one. It is shown that taking into 
acount non-stationarity of radiative transfer in Ly-a line can lead to an apreciable additional 
acceleration of primordial hydrogen recombination. 
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